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£NJ ■ Abstract 

In this paper we present a methodology employing statistical analysis and stochastic geometry to 
study geometric routing schemes in wireless ad-hoc networks. In particular, we analyze the network layer 
performance of one such scheme, the random 1/2-disk routing scheme, which is a localized geometric 
routing scheme in which each node chooses the next relay randomly among the nodes within its 
transmission range and in the general direction of the destination. The techniques developed in this paper 
enable us to establish the asymptotic connectivity and the convergence results for the mean and variance 
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■<^j- ■ of the routing path lengths generated by geometric routing schemes in random wireless networks. In 
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particular, we determine the sufficient conditions that ensure the asymptotic connectivity for both dense 
and large-scale ad-hoc networks deploying the random 1/2-disk routing scheme. Furthermore, we show 
that the expected length of the path generated by the random 1 /2-disk routing scheme normalized by the 
length of the path generated by the ideal direct-line routing, converges to 37r/4 asymptotically. Moreover, 
we show that the variance of the routing path length normalized by its expected value converges to 
97r 2 /64— 1 asymptotically; this indicates that the dispersion of the individual routing-path lengths around 
their mean remains constant relative to their mean regardless of the granularity and size of the network. 
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I. Introduction 

A wireless ad-hoc network consists of autonomous wireless nodes that collaborate on com- 
municating information in the absence of a fixed infrastructure. Each of the nodes might act as a 
source/destination node or as a relay. Communication occurs between a source-destination pair 
through a single-hop transmission if they are close enough, or through multi-hop transmissions 
over intermediate relaying nodes if they are far apart. The selection of relaying nodes along the 
multi-hop path is governed by the adopted routing scheme. 

The conventional method to establish a routing path between a given source-destination pair 
is through exchanges of control packets containing the complete network topology information 
JT), which creates scalability issues when the network size becomes large. One way to reduce 
the overhead for global topology inquiries is to build routes on demand via flooding techniques 
0. However, such routing protocols essentially suffer from a similar issue of large signaling 
overheads. To deal with the above issues, Takagi and Kleinrock O introduced the first geo- 
graphical (or position-based) routing scheme, coined as Most Forward within Radius (MFR), 
based on the notion of progress[]: Given a transmitting node S and a destination node Dst, 
the progress at relay node V is defined as the projection of the line segment SV onto the line 
connecting S and Dst. In MFR, each node forwards the packet to the neighbor with the largest 
progress (e.g., node V2 in Fig. [T|), or discards the packet if none of its neighbors are closer to 
the destination than itself. There are some other variants of the geographical routing scheme in 
the literature Il4~lll5l lfl2l . which are similar to MFR. In [4], the authors introduced the Nearest 
Forward Progress (NFP) method that selects the nearest neighbor of the transmitter with forward 
(positive) progress (e.g., node V\ in Fig. [D; in 0, the Compass Routing (also referred to as the 
DIR method) was proposed, where the neighbor closest to the line connecting the sender and 
the destination is chosen (e.g., node V 3 in Fig. [[]); in [fl2|. the authors considered the Shortest 
Remaining Distance (SRD) method, where the neighbor closet to the destination is selected as 
the relay (e.g., node V4 in Fig. [T]). 

Geographical routing protocols might fail for some network configurations due to dead-ends 
or routing loops. In these cases, alternative routing strategies, such as route discovery based 

'it should be noted that the reduction in complexity comes at the cost of knowing the location of the neighboring nodes in 
addition to that of the destination. 
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Fig. 1. Some variants of geometric routing schemes: The source node S has different choices to find a relay node for further 
forwarding a message to the destination Dst. V\ = Nearest Forward Progress (NFP), V2 = Most Forward within Radius (MFR), 
Vz = Compass Routing (DIR), V4 = Shortest Remaining Distance (SRD). 

on flooding [7] and face routing (§1 can be deployed. However, it has been shown in [9] that 
for dense wireless networks, the MFR-like routing strategies will succeed with high probability 
and there is no need to resort to recovery methods such as face routing. In this paper we study 
the network layer performance of geographical routing schemes in such dense or large wireless 
networks; and we expect to observe a similar high-probability successful routing performance 
(the proof of this claim is presented in Section IIV-BD . 

Below we present a methodology employing statistical analysis and stochastic geometry to 
study geometric routing schemes in wireless ad-hoc networks. We consider a wireless ad-hoc 
network consisting of wireless nodes that are distributed according to a Poisson point process over 
a circular area, where nodes are randomly grouped in source-destination pairs and can establish 
direct communication links with other nodes that are within a certain range. We determine the 
conditions under which, in such a network, all source-destination node pairs are connected via the 
adopted geographical routing scheme with high probability and quantify the asymptotic statistics 
(mean and variance) for the length of the generated routing paths. In particular, we focus on a 
variant of the geographical routing schemes, namely the random 1/2-disk routing scheme, as 
an example, where each node chooses the next relay uniformly at random among the nodes 
in its transmission range over a 1/2-disk with radius R oriented towards the destination. This 
scheme is similar to the geometric routing scheme discussed in Q, in which one of the nodes 
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with forward progress is chosen as a relay at random, arguing that there is a trade-off between 
progress and transmission success. 

We chose the random 1/2-disk routing scheme mainly for tractability and simplicity in 
mathematical characterization. However, the solution techniques developed in this paper can 
be used (with some modifications) to study other variants of geographical routing schemes, 
such as MFR, NFP, DIR, etc, which will be further discussed in Section |VIJ Moreover, the 
random 1/2-disk routing scheme can be used to model situations where nodes have partial or 
imprecise routing information and the locally optimal selection criterion of greedy forwarding 
schemes fails [6], e.g., when nodes have perfect knowledge about their destination locations but 
imprecise information about their own locations, or when nodes only know the half-plane over 
which the final destination lies such that randomly forwarding the packet to a node in the general 
direction of the destination is a plausible choice. 

There has been a considerable interest regarding the network connectivity and the average 
length of the route generated by geographical routing schemes under different network set- 
tings ll6ll [fT0ll -[[T4 l| . The authors in ifTOl considered a wireless network that consists of n nodes 
uniformly distributed over a disc of unit area with each node transmission covering an area 
of r(n) = (log n + c(n))/n. They show that this network is connected asymptotically with 
probability one if and only if c(n) — > oo as n — > oo. Although the asymptotic expression that 
they derived for the sufficient transmission range is similar to ours, their notion of connectivity 
is quite different from ours. In ifTOll . the network is connected as long as it is percolated, i.e., the 
network contains an infinite-order component, where no constraints are considered for the paths 
connecting source-destination pairs. However, the routing paths that we consider in this work have 
more structure such that we need a different proof technique to prove the asymptotic connectivity 
of the network. Xing et al. showed in ifTTI that the route establishment can be guaranteed 
between any source-destination pair using greedy forwarding schemes if the transmission radius 
is larger than twice the sensing radius in a fully covered homogeneous wireless sensor network. 
In lfT2l the authors derived the critical transmission radius to be \J ' l3ol ° gn which ensures network 
connectivity asymptotically almost surely (a.a.s.) based on the SRD routing method, where 
Po = 1/(2tt/3 - y/3/2). 

In lfT3ll . Bordenave considered the maximal progress navigation for small world networks 



5 

and showed that small world navigation^ is regenerative. It is shown furthermore in [13 ] that 
as the cardinality of the navigation (or routing) path grows, the expected number of hops 
converges, without providing an explicit value for the limit. Baccelli et al. IfPfll introduced a time- 
space opportunistic routing scheme for wireless ad-hoc networks which utilizes a self- selection 
procedure at the receivers. They show through simulations that such opportunistic schemes 
can significantly outperform traditional routing schemes when properly optimized. Furthermore, 
they analytically proved the asymptotic convergence of such schemes. In J6l, Subramanian and 
Shakkottai studied the routing delay (measured by the expected length of the routing path) 
of geographic routing schemes when the information available to network nodes is limited or 
imprecise. They showed that one can still achieve the same delay scaling even with limited 
information. Note that the asymptotic delay expression derived in [0 is similar to the one 
we derive in this paper; however, our proof technique is more constructive and enables us to 
derive tight bounds for the mean and the variance of the routing-path lengths in a network of 
arbitrary size, together with the exact expressions for their asymptotes. Moreover, in [|6|| the 
authors assume a continuum model for the sensor network and presumes that the progress (as 
defined in and described earlier) at nodes along the routing path form a sequence of i.i.d. 
random variables. However, as we show later (cf. Proposition [D, this assumption may not hold 
for Poisson distributed networks of arbitrary finite sizes as the distribution of nodes contained 
in the transmission range of a given node along a routing path depends on the history of the 
routing path up to this node, i.e., the progress at each hop is history dependent. Hence, it is 
neither independent nor identically distributed; but, as the size of the network (either density or 
area) goes to infinity, the distribution of the sequence of progresses along the routing path, in 
fact, converges to an i.i.d. sequence of random variables. 

The remainder of this paper is organized as follows. In Section [II] we introduce the system 
model and describe the random 1/2-disk routing scheme. Then we define the notion of con- 
nectivity based on generic geometric routing schemes and state the main results of the paper 
in a theorem regarding the connectivity and the statistical performance of the random 1/2- 
disk routing scheme. In Sections [Till and [IV] we prove the claims made in this theorem. In 
Section [Till we establish sufficient conditions on the transmission range that ensure the existence 



2 This routing scheme, unlike ours, assumes nonnegative progress in each hop. 
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of a relaying node in every direction of a transmitting node for both dense and large-scale 
networks. In Section [TV] we study the stochastic properties of the paths generated by the 
random 1/2-disk routing scheme. Specifically, in Section |IV-A[ we show that the process of 
path establishment by the random 1/2-disk routing scheme can be approximated by a Markov 
process that converges (statistically) to the actual process asymptotically. In Section ITV-B I using 
the Markov characterization, we derive the asymptotic expression for the expected length, and in 
Section IIV-CI we derive the asymptotic expression for the variance of the length of the random 
1/2-disk routing paths. In Section|Vl we present some simulation results to validate our analytical 
results. In Section [VH we present some guidelines on how to generalize the results derived for 
the random 1/2-disk routing scheme to other variants of the geometric routing schemes. We 
conclude the paper in Section IVIIl 

II. System Model 

Consider a circular area A over which a network of wireless nodes resides^. Nodes are 
distributed according to a homogeneous Poisson point process with density A. Each node picks 
a destination node uniformly at random among all other nodes in the network, and operates with 
a fixed transmissionpower that can cover a disk of radius R = R(X, \A\), where \A\ denotes 
the area of region AZ- 

For a generic geometric routing scheme, when the targeted destination node is out of the one- 
hop transmission range (R) of a given transmitting node, the next relay is selected (based on 
some rules) among the nodes contained in the relay selection region (RSR) of the transmitting 
node, where the RSR, in general, can be any subset of a full disk of radius R centered at the 
transmitting node. For example, the RSR for all the geometric routing schemes cited in the 
introduction section is a 1/2-disk of radius R centered at the transmitting node and oriented 
towards the destination (denoted by |RSR). We define the rule that governs the selection of 
the next relay in each node's RSR as the relay selection rule (RSL). For example, the RSL for 
MFR is to choose the node with the largest "progress" towards the destination among the nodes 

3 The results will carry over, with some minor considerations, to any convex region with bounded curvature. 

4 As mentioned earlier, we are only interested in the network layer performance of the network; as such, we do not consider 
physical layer related issues such as interference. However, as a rule of thumb (cf. ||9)), to minimize the interference among 
wireless nodes we are interested in the smallest transmission radius that ensures network connectivity in this paper. 
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contained in its |RSR. We define the progress x' v at a relay node V as in [3], and described in 
the introduction section. 

We define the network to be connected if for any source-destination node pair in the network, 
there exists a path constructed by & finite sequence of relay nodes complying with the RSL, with 
high probabilitj^; henceforth, we call such a relay sequence a routing path. Note that a node can 
potentially act as a relay only if it is contained in the RSR of the current transmitting node. For 
the sake of definition, we claim that the network is connected if the network node set is empty. 
In this paper we study a special case of localized geometric routing schemes, namely the random 
1/2-disk routing scheme, where for each transmitting node S in the network, as illustrated in Fig. 
|2l the next relay V is selected uniformly at random among the nodes contained in the |RSR of 
S. We denote the relay selection rule of the random 1/2-disk routing scheme by rRSL. Observe 
that according to our routing scheme, the next chosen relay might be farther away from the 
destination than the current transmitting node. 

In the following, we present a theorem that summarizes the main result of this paper on 
the random 1/2-disk routing scheme, regarding i) the sufficient conditions on R(X, \A\), which 
ensure the existence of a relaying node in any direction of a particular transmitting node based 
on a generalized version of |RSR; ii) the mean asymptotes of the path-lengths established by 
the random 1/2-disk routing scheme; iii) the corresponding variance asymptotes; and iv) the 
asymptotic network connectivity with the random 1/2-disk routing scheme. For the generalized 
version of the |RSR, we assume that the RSR of a node is a wedge of angle 27r?7 with radius 
R, where < rj < 1 (hereafter called //-disk or 7/RSR, interchangeably). Hence, the |RSR is a 
special case of the 77RSR with 77 = 1/2. 

For notational convenience, we let N := A|A| designate the expected number of nodes in the 
network region of area \ A\ and d = d(N) := denote the normalized area of a full disk with 
radius R relative to the area of the whole region, such that dN is the expected number of nodes 
in such a disk. The asymptotic nature of the results presented in this paper is due to N — > 00, 
which can represent results for either large-scale networks (i.e., when \A\ — > 00 with a fixed A) 
or dense networks (i.e., when A — > 00 with a fixed \A\). Also, f(n) = O (g(n)) means that there 
exist positive constants c\ and M such that f(n)/g(n) < c\ whenever n > M, f(n) = o(g(n)) 



5 According to this definition, the network is connected if starting from any source and choosing relays based on the routing 
scheme, the destination is reachable with high probability. 
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Fig. 2. The random 1/2-disk routing scheme. 

means that lim f(n)/g(n) — > as n — > oo, f(n) ~ g(n) means that lim f(n)/g(n) — > 1 as 
n — >■ oo, and /(n) = (g(n)) means that there exist positive constants c\, c 2 and M such that 

ci < f(n)/g(n) < c 2 whenever n > M. 

Theorem 1. Consider a Poisson distributed wireless network with an average node population 
N deployed over a circular area A. Assume all nodes have the same transmission range R{N) 
that covers a normalized area d = d(N) and let x' be the progress at each node. Then: 

i) the rj-disk of each node in the network pointing at any direction in which its targeted 
destinations may lie contains at least one relaying node a.a.s., if d = o (iV^ 2 / 3 ) and 
rjdN + log d — > +oo as N — > oo; 

ii) given (i), the length v of the random 1/2-disk routing path is a.a.s. finite with the asymptotic 
expected value converging to ff^ as ^ ~^ 00 » specifically, the expected length of the 
random 1/2-disk routing path connecting a source-destination pair that is h-distance apart 
satisfies E{v \ h) ~ h/E(x') = as N -> oo; 

iii) the variance of the routing path length normalized by its mean, Variv) /E[y), converges 



iv) given (i) {and consequently (ii)J, the network is connected a.a.s. with the random 1/2-disk 
routing scheme. 




Proof: Here we only sketch the outline of the proof and present the respective details in 



the following sections. In Section [TTT], we show that for random networks, choosing R(N) such 
that d = o (iV -2 / 3 ) and rjdN + log d — > +00 as N — > 00 guarantees the existence of at least 
one relaying node in the //-disk of each network node pointing at any directions in which their 
targeted destinations may lid_| a.a.s.. To this end, we first derive an upper bound cr(iV) on the 
probability that the //-disk of some nodes in the network pointing at some directions is empty. 
Then we show that choosing d(N) as mentioned before ensures the asymptotic convergence of 
a(N) to zero as N — > 00. This ensures the existence of a relaying node in every direction of 
a particular transmitting node and ascertains the possibility of packet delivery to a particular 
destination from any direction. 

In Section [TV] given the existence of a relaying node in every direction of a particular 
transmitting node, we show that the length of the random 1/2-disk routing path connecting 
a source to its destination (that are /i-distance apart) is finite almost surely. This shows that 
starting from a source and following the random 1/2-disk routing scheme we can reach the 
destination in finitely many hops a.a.s. (regardless of the specific realization for the network 
or the routing path); hence the network is connected with the random 1/2-disk routing scheme 
a.a.s.. More specifically, we show that (in Proposition [I]), we can approximate the process of 
construction or formation of a routing path between a source-destination pair as a Markov process 
that converges to the actual process asymptotically. Using this characterization, we then derive 
the asymptotic expressions for the mean and variance of the routing path length generated by the 
random 1/2-disk routing scheme between a source-destination pair that is /i-distance apart and 
show that they are asymptotic to = j| and ^ r jf^ E [y) — — lj E (u), respectively. 



III. Theorem [Hi Proof: Uniform Relaying Capability 

In this section we derive the sufficient conditions on R(N) that ensures, for any node in 
the network, its 7/-disk pointing at any directions over which its targeted destinations may lie 
contains at least one potential relaying node. To this end, we first characterize the upper bound 
on the probability a(N) that, for some network nodes, there are certain directions at which their 
//-disks are empty; we then choose R such that this bound is vanishingly small. In this process, 

6 A specific node might act as a relay for multiple source-destination pairs. 
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we can distinguish between two types of network nodes based on their distances to the edge of 
the network: Nodes that are farther than R away from the edge of the network, which we call 
interior nodes, and nodes that are closer than R to the edge of the network, which we call edge 
nodes. 

For interior nodes, it is clear that the node distribution in their ry-disks, pointing at any direction, 
is the same. Therefore, the existence probability of an empty 7/-disk for an interior node is 
independent of its targeted destination direction. However, due to the proximity of edge nodes 
to the boundary of the network, the existence probability of an empty 7/-disk for an edge node 
highly depends on its destination orientation. For example, the 77-disks that fall partly outside 
the network region are more likely to be empty than the ones that are fully contained in the 
network region. Hence, we derive the probabilities of a node having an empty //-disk in some 
direction separately for the interior nodes and the edge nodes, denoted by cr'(N) and a"(N), 
respectively. 

Recall that a 7/-disk is a wedge of angle 2Tn] and radius R, with < 77 < 1. Hence, the 
TjRSR is a special case of 7/-disk with 77 = 1/2. Each //-disk has an expected number of nodes 
r]dN. As shown in Section IIII-CL the existence probability of an empty 77-disk increases as 77 
decreases. However, we can show that the expected length of the routing path connecting a 
source to its destination will decrease as 77 decreases. Hence, there exist a tradeoff between the 
existence probability of an empty 77-disk (i.e., a disconnected node) and the expected length of 
the routing path between a source-destination pair parameterized by r\. We leave the study of 
this tradeoff to a future work and only derive (in Section ITVl) the mean and variance of the path 
length connecting a source-destination pair when 77 = 1/2. 

A. Calculation of o~'(N) 

Consider an interior node x, fixed for now. Given i > 1 nodes are in the transmission range 
of x, their directions in reference to x are independent and uniformly distributed on [0, 2tt]. The 
probability that x has an empty 77-disk in some direction equals the probability U^rj) that the 
angle of the widest wedge containing none of these i nodes is at least 27771". It is not difficult to 
give a simple upper bound on Ui(rj): Of the i nodes, without loss of generality (W.L.O.G.), we 
can assume that (at least) one is at one edge of an empty wedge with angle of 27771% while the 
other i — 1 are distributed independently and uniformly in the remainder of the full transmission 
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disk, as shown in Fig. [3j Hence, we obtain Ui{rj) < z(l — r]) % l , for i > 1. Of course, if i — 
the probability is t/ (?7) = 1- 




Fig. 3. A realization for which the widest wedge between the nodes is of an angle at least 2r]-K. 



One can obtain a more precise expression for Ui{rj) using results in 031, page 188: 

min{LV?7j,i} , 

for i > 1, where [aj is the largest integer smaller than a. This expression is based on the 
inclusion-exclusion principle for the probability of the union of events, for which the first term 
in the sum provides an upper bound and the first two terms provide a lower bound. 

Let L := yJ\A\/it = R/Vdbe the network radius and a\ := ir(L — R) 2 /\A\ be the normalized 
area of the network interior. Then, averaging over i (number of the nodes in the transmission 
range of x) and over the number of interior nodes, we have: 

k — 1 j — 1 i — 

1 rlN 2 p- r l dN f 1 \ 

where the first inequality is due to union bound and the second inequality is due to the fact that 

(1 - (1 - a^Y 1 < a^ 1 for k > 1. 
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B. Calculation of a" (N) 

So far we have considered the interior nodes that are at least i?-distance away from the 
boundary of the network region. Now, we consider edge nodes that are within R of the network 
edge. Therefore, some ??-disks of an edge node may fall partially (up to half) outside the region, 
which increases the chance that they are empty. We refer to this phenomenon as the edge effect. 
Since the network region is circular, the number of such edge nodes equals (2 — yfd)\fdN, which 
is of order 6 ys/dN^j. We need to determine how their contribution to cr(N) differs from the 
interior nodes. 




Fig. 4. Edge curvature. 

Consider an edge node e, (5'i?)-distance away from the network edge, with < 8' < 1. 
As shown in Fig. HI we take node e as the pole and the ray ev (perpendicular to the network 
edge) as the polar axis of the local (polar) coordinates at node e. We observe that, due to the 
curvature of the network edge, the overlap of node e's transmission range with the network 
region is larger than what it would be if the network's edge were straight (i.e., the line passing 
through the intersection points A and B in Fig. HI). This area-difference (the shaded area in Fig. 
B} is no larger than \fdR 2 containing an expected number of nodes on the order of (<i 3,/2 ./V), 
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where the maximum area-difference is obtained when node e is located on the straight network- 
edge approximation line (i.e., in the middle of AB in Fig. H]). Accordingly, we make a further 
simplifying assumption that d = o (iV -2 / 3 ) ; this is equivalent to a practical assumption that the 
ratio between the transmission range R and the radius L of the network region goes to zero 
fast enough such that the expected number of nodes in the shaded area of Fig. @] goes to zero 
as N — y oo. Then the error in calculating the probability of any event in the following will 
be a factor of no more than e kd,i/2N — > 1, where A; is a finite constant Henceforth, for large 
N, we proceed as if the network region is straight wherever it intersects with an edge node's 
transmission disk, i.e., we neglect the effect of such shaded areas as shown in Fig. HI 



Network Edge 




Fig. 5. Intersection of the 77-disk with the network region. 

We argued in the beginning of this section that, for edge node e, the probability of an 77-disk 
being empty, depends highly on its orientation. Let us consider this claim more closely. Let 
if = cos -1 (5) G (0, 7r/2), as shown in Fig. [5j where 5R is the distance between node e and the 
straight approximation of the network edge as defined before with < 6 < 1. Note that all the 
77-disks are oriented towards the destination node. Hence, for all ?7-disks that are oriented towards 

7 This rate will apply as long as the region is convex with finite/smooth curvature. 
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an angle in the range (—ip, f), we must have that the destination is within node e's transmission 
range. Therefore, we only need to be concerned with empty 77-disks oriented towards an angle in 
the range (ip, 2% — ip). The //-disks oriented to an angle in the range (—p — 7]Tr, —p) U (ip, ip + rjn) 
are partially outside the network region, as illustrated in Fig. (5J and those oriented to any angle 
in (ip + r/7r, 27r — ip — r]ir) are fully contained inside the network region. Note that here, all the 
angles are measured relative to the polar axis ev . 

We now compute cr"(N) for node e. Let a 3 := n(L 2 - (L - R) 2 )/\A\ = \fd(2 - \fd) and 
a 2 := ir(L 2 — (L — 2R) 2 )/\A\ = A\/d(l - y/d) be the normalized areas of the network edge 
region and the network extended edge regioro respectively. First, suppose that there are no nodes 
within the transmission range of node e; this event occurs with probability no greater than 

1 v^-o^MO' 



ti « l-(l-S)'feWV 11 J J(i 2a 2 J 



-a 2 AT 

Z=l * V " a 2 y J 

, 2a 3 iVe- dJV / 2 _ 2(2 - ^d)VdN e - dN / 2 

— I — e -a 2 N ~ j _ e -4(l- v / d)rf 3 / 2 A r ' ^ 

where again the inequalities are due to the union bound and the fact that a 3 /a 2 > 1/2 and 
a 3 /a 2 -»■ 1/2 as iV ->• 00, such that (1 - (1 - asM)')" 1 < C 1 - 2 ^) _1 < 2 - 

Second, suppose that there are i > 1 nodes in the intersection of node e's transmission range 
with the network region. If an empty 77-disk exists and it is completely contained within the 
network region, W.L.O.G., there should be a node on its left edge at some angle 9 G (cp + 
2t]ti, 2-K — ip). However, for an empty 77-disk that is partially contained within the network region 
there should be, again W.L.O.G., a node at an angle 9 E (tp + rjir, (p + 2r)n) or 9 £ (—ip, —Lp + r/Tr) 
on the left edge of the 77-disk (note that, as discussed earlier, no 77-disks can be oriented towards 
an angle in (—ip, ip)). Clearly, the existence probability of such empty 77-disks (that is partially 
contained in the region A) increases as either 8 or \9\ decreases. The area of the intersection 
between such an 77-disk (that is partially contained in the region A) and the network region A is 
that of a wedge with angle \9\—<p (wedge AeB in Fig. [5j plus a triangle abutting the right edge 
of the wedge (triangle BeC in Fig. [5j. In fact for an arbitrary small e, if either 8 > sin(3e7r) or 
9 > ip + ?77r + 2e7r, the area of the intersection between the 77-disk and the network region is at 
least (77/2 + e)TiR 2 . Otherwise, it is at least r)nR 2 /2. Thus, averaging over 8, 9 and the number 

8 The extended edge region is the area of the network that is within 2R of the network edge. 
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of edge nodes, the probability that some edge nodes have empty 77-disks in some directions, 
a"(N), is derived to be no more than 
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7/ + 2e ^- 1 + (1-2^(1 



+ 2^(1 



2// 



l + 8e' 



l + 8e' 



+ 



< 



2r/(l-(r / /2 + e)) l " 1 + (l-2r/)(l-r / ) 



2d 3 / 2 iV 2 



;i- e -^)(i-^ 



I2ixe e i+s^ + 67ree i+s<= + 37ree 1+^ + 2e 2 _|_ e ^ aJV 



(3) 



for arbitrary e > 0. Choosing e = 2lo ^ N ' , together with ©, yields the upper bound for the 
probability that some edge nodes has an empty r/-disk oriented in some direction: 



a" (iV) < 



400vr (log dKf 
Vd 



n dN 16(dN) 



-e^ dN + 4VdNe^ dN , 

y/d 



(4) 



for large enough dN where the last summand is the probability that some edge nodes have no 
other nodes within their transmission ranges, derived in ©. 



C. Calculation of a(N) 

Finally, summing (OQ) and ©, we obtain the bound cr(N) on the probability that some nodes 
in the network have empty 77-disks looking in some directions as: 

a(N) < 40071 ^ dN )\~W + 1Q (dN)\. vdN + ^Ne*™ + AdN 2 e~ vdN . (5) 
yd yd 

This bound on er(iV) is asymptotic to 4007r (^s dAr ) e ~2 dN , which goes to zero if rjdN+logd —> 
00 as 00. Hence, setting d = clo % N with c > I/77, we obtain that every node in the network 
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have at least one relaying node in every direction over which their targeted destinations may lie 
with probability approaching one as N — > oo, which shows the consistency between our result 
and the ones derived in IfTOl , |[T6l and IfTTl for r] — 1. 

Remark 1. Setting d = is equivalent to setting R(X, \A\) = logA+ A log|A| /or c > 1/rj. 

In particular, for the case of dense networks (i.e., X — > oo with a finite \A\) and for the case of 
large-scale networks (i.e., \A\ — > oo with a finite X), setting R(X) = K^/\ogX/X and R(\A\) = 
K y^log \A\ respectively, with a large enough constant K, guarantees the existence of relaying 
nodes in a "uniform" manner around each node in the network. 

IV. Theorem \\\ii-iv Proof: Path Length Statistics and Connectivity 

Assuming that each network node has at least one relaying node in every direction, we now 
investigate the question of how long the path generated by the random 77-disk routing scheme 
is, where we focus on the 77 = 1/2 case in this paper. To answer this question, we need to 
characterize the process of path establishment (from a given source to its destination) by the 
random 1/2-disk routing scheme. 

Consider an arbitrary source-destination pair that is /i-distance apart. We set the destination 
node at the origin and assume that the routing path starts from the source node at X = (—h, 0), 
where X n is the (Cartesian) coordinate of the n th relay node along the routing path and r n := 
\\X n \\ is the (Euclidean) distance of the n th relay node from the destination. 

More specifically, the routing path starts at the source node X = (—h,0) with its |RSR -Do 
that is a 1/2-disk with radius R centered at X and oriented towards the destination at (0,0). The 
next relay X\ is selected at random from those contained in D (the rRSL rule). This induces a 
new iRSR Di, also a 1/2-disk but centered at X± and oriented towards the destination. Relay X 2 
is selected randomly among the nodes in D\, and the process continues in the same manner until 
the destination is within the transmission range. We claim that the routing path has converged 
(or is established) whenever it enters the transmission/reception range of the final destination, 
i.e., r T < R, for some r G {1, 2, • • ■ }. In Fig. [6l we illustrate the progress of routing towards 
the destination. 

Define the routing increment as F n+1 := ||X n || — ||X n+1 || = r n — r n+1 , and let (j)(D n ) be the 
number of nodes in D n . In the next section we investigate how similar Y n (and consequently 
r n ) is to a Markov process. 
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A. Markov Approximation 

Recall the definition of the routing increment Y n . Note that even though the underlying 
distribution of the network nodes is Poisson and the new relays are chosen uniformly at random 
within each |RSR, the increments Yi, Y 2 , ■ ■ ■ are neither independent nor identically distributed. 
This is due to the fact that the orientations of all |RSRs are pointing to a common node 
(destination) and might overlap, as shown in Fig. [6l More specifically, the overlap of D n with 
some Dj, < j < n, results in the dependence of the spatial distribution of nodes in D n (and 
consequently Y n+ i), not only on X n , but also possibly on Xj, < j < n. This dependence 
increases as the packet gets closer to the destinationj. In addition, due to the overlap of D n with 
D n -\ (and perhaps some Dj, < j < n — 1), the nodes contained in D n are not uniformly 
distributed over D n as one would expect for a Poisson distributed network (cf. Proposition [Q). 
As such, the process of path establishment (from a source to its destination) by the random 
1/2-disk routing scheme is not a Markov process; however, as shown in Proposition [fl it can 
be approximated by a Markov process that converges (stochastically) to the actual process as 
iV ->■ oo. 

Since tracking the dependence of Y n+ i on all Xj, j < n, is extremely tedious, in this 

'Because the overlapping area between D n and D n -i, D„_2, ■ • • increases as the packet gets closer to the destination 
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work we only investigate the dependence of Y n+1 on X n and X n _ x , and neglect the effect 
of X n _ 2 , X n _ 3 , ■ ■ ■ which is justifiable under certain conditiono- In particular, we determine 
how close the distribution of the nodes in D n is to a uniform distribution when just knowing 
the locations of the current and previous nodes X n , X n _i, this is equivalent to showing how 
similar {r n } is to a Markov process of order one. 

Note that conditioned on X n (or equivalently on D n ) and the existence of a relaying node in 
every direction of X n , <p(D n ) is Poisson distributed with intensity A|-D n | and zero mass at zero. 
What is less clear, however, is the nature of (f>(D n ), given X n _ x or equivalently D n _\. Also 
note that throughout this section we assume that the |RSR of each network node looking in 
any direction is nonempty, i.e., 4>{D n ) > 0, for all n and all source-destination pairs and every 
path between them, unless otherwise stated. We emphasize that, in what follows, conditioning 
on 4>(D n ) > means we only know that there is at least one node in D n ; however, conditioning 
on 4>(D n ) means we know the exact number of nodes in D n . 

Observe that D n only depends on X n . Given X n , X n _i, <j)(D n _i), and 4>{D n ) > 0, the number 
of nodes in D n _ x D n := D n _i n D n is 0(iViA») ~ Binomial - 1, '^^i"' ) + 

1{x„_i6D»} an d independent of the number of nodes in D^_ x D n , which is 0(D^_ 1 D n ) ~ 
Poisson(A|D^_ 1 D n |), where C c := A — C denotes the complement of C with respect to 
network region A and 1{.} represents the indicator function, i.e., = 1 if the event in 
the subscript happens and = otherwise. Moreover, conditioned additionally on the two 
random variables <p(D n _iD n ) and 0(-D^_ 1 Z) n ), each collection of nodes (located in D n _iD n 
and D^^Dn) is uniformly distributed on the respective areas. This does not, however, imply 
that the combined collection of nodes is uniformly distributed on D n . The combined points are 
uniformly distributed on D n only if the (conditional) expected proportion of points in D n _iD n 
is E ( % D ^ n) | 0(AO > 0,^(iVi) > 0,X n ,X n _!) = 

Nonetheless, according to the following proposition, the error resulted from proceeding as 
if X n+ i is located uniformly on D n is negligible for large N . Essentially, knowing X n , the 
distribution of nodes in D n is almost uniform over D n and independent of the location of the 

10 The analysis gets more complicated as we consider a longer history of the previous relaying nodes that their RSRs intersect 
with D n , i.e., X n -2, X n s, ■ ■ ■ ,X n _fc, but it can be shown that if k = o yy/dN^j, the error resulting from neglecting the 
previous relaying nodes should remain in the order of O (l/(dN)) where l/(dN) is the error resulting from neglecting X n -i, 
as shown in Proposition Q] 
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previous relaying node X n _i for large N. 

Proposition 1. Assume that every node in the network has at least one relaying node in all 
directions^ and the locations of current and previous relay nodes, X n and X n _i, are given. Then 
the distribution of the nodes located in D n (the ^RSR of the current node) converges to a uniform 
distribution as N — )■ oo. In particular, the conditional probability of selecting the next node 
X n+l from D n ^D n , i.e., p(X n _ u X n ) := E ( % D ^ n) | <j>{D n ) > 0, 0(Z> n _i) > 0, X n , X n ^ 
satisfies 

2 / x -a 2 (n)dN\ \D n -lD n \ _ fv ^ , ^ \D n _xD n \ 



where ai(n) > 2 and < ^(n) < 1 are independent of N. 

Proof: Refer to Appendix |Al ■ 
Observe that according to ©, given the location of the previous relay node X n -i, it is less 
likely that the next relay X n+ i is selected from D n ~\D n as opposed to the case where the nodes 
were actually uniformly distributed in D n . However, we have p(X n -i,X n ) — > \D n D n _i\/\D n \ 
as N — > oo. Hence, we obtain that for large N, Y n+ i only depends on X n and is (almost) 
independent of X n -\. In other words, the routing increment Y at the current relay is only a 
function of the current node and is independent of the history of the routing path for large N . 
Nevertheless, Yi,Y 2 , ■ ■ ■ are not identically distributed and as shown in the next section, Y n+ i 
is in fact a function of r n . As such, for large N we can proceed as if the process that governs 
the path establishment by the random 1/2-disk routing scheme is a non-homogeneous Markov 
process^. 

B. Theorem and iv Proof: Expected Length of the Random 1/2-Disk Routing Path and 
Network Connectivity 

According to Section IIV-AI we can approximately model the distance evolution {r n } of the 
routing path from a source node to its destination node as a Markov process solely characterized 

"Note that by TheoremQ] the sufficient condition for this to happen is rjdN + log d — > +00 as N — > 00, which implies that 
dN — > 00 and d — > for smallest transmission radius ||9)- 

12 In this section and what follows, we ignore the edge effect. More precisely, we assume that D n n A ~ D n irrespective of 
the location of X n . 
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Fig. 7. Distance between the next relay and the current node projected onto to the local coordinates at the current node. 

by its step sizes {Y n }. Let (x' n+1 ,y' n+1 ) be the projection of X n+1 — X n onto the local Cartesian 
coordinates with node X n as the origin and the z-axis pointing from X n to the destination node 
as shown in Fig. [7J Hence, 

r n+l = \j Vn ~ X'n+l) 2 + Vn+1 , ( 7 ) 

characterizes the distance evolution of the routing path at the n th hop. According to Proposition 
HI X n+ i is uniformly distributed on D n for large enough N; hence {(x' n , y' n )} is an i.i.d. sequence 
of random variables with < x' n < R and — R < y' n < R for all n, whenever N is large enough. 

Define Ur '■= inf{n : r n < r, r = h}, r > R, to be the index of the first relay node (along 
the routing path) that gets closer than r to the destination when the source and destination nodes 
are /i-distance apart. Hence, represents the first time the routing path enters the reception 
range of the destination and + 1 quantifies the length of the routing path. It is easy to show 
that vi^ is a stopping time lfP9l and 

r — R < r (h) < r . 

Furthermore, let g(r, x' , y') := \/Jr~— x') 2 + y' 2 — r. Observe that g is a nonincreasing function 
over r > R, for fixed (x', y'), and g(r n , x' n+1 , y' n+1 ) = —Y n+ i. Thus, for n < uf 1 ', we have r„ > r 
and 

—x n+l < r n+ i — r n = g(r n , x n+1 , y n+1 ) < g(r, x n+1 , y n+ x) ■ 
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Hence, for a source-destination pair that is /i-distance apart (r = h), we have 

r- R< r v ( h ) <h + ^g(r, x' n , y' n ) , (8a) 

n=l 

h + J2(- X 'n)^ r 4 h) - V - (8b) 

n=l 

Note, as well, that (refer to Appendix |B]) 

~ = E (-<) <E («?(r, a4, j/J) < E (<?(#, x' n , j/J) < < . (9) 

07T 4 



Now, applying Wald's equality 112011 to (l8al) and (l8bl) and rearranging, we obtain a bound on 
the expected value of the stopping time Ur . 

3 ^ h ~ r ) < E f v W I M < fe-r + i? 4^ 

4i2 — I / — _ E (^(r,x' n ,^)) " -E(^(r,x' n ,<)) " J? " 

Substituting r with i? we obtain a general bound for the expected length of routing path (minus 
one) between a source-destination pair that is h distance apart as 

T(H SE (* ,|k ) S TT 

This implies that the length of the random 1/2-disk routing path is almost surely (a.s.) finite when 
each network node has at least one node in its |RSR looking in any direction, which happens with 
probability no less than 1 — cr(N) as obtained in Q. In other words, when dN/2 + logrf — > oo 
as N — > oo, we obtain that Pr (v^ < ooj — > 1 as N — > oo. This in turn shows that given 
dN/2 + logd — > oo as N — > oo, every path starting from any source will reach its destination 
in finitely many hops a.a.s., which proves that the network is connected employing the random 
1/2-disk routing scheme, according to the connectivity definition in Section HH 

When the ratio h/R (i.e., the ratio between the source-destination distance and the transmission 
range) is large, we can obtain a tighter bound on the expected length of the routing path between a 
source-destination pair with h separation. For the following, we assume h > 2R. Since r^ h ) < r, 
we must have E \ hj < E \ hj +E (v^ \ rj . Thus, by (flOl) and proper substitutions, 



we have 



E (v$> | h\ 



h/R ~ -E(g(r,x' n ,y' n )) h 



■ 
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for all R < r < h. Using 



-x' n <g(r,x' n ,y' n ) < -x' n + 



I \2 



2(r - R) 



(ID 



and d23B we get E (g(r, x' n , y' n )) < -ff + 



R 2 



8{r-R) 
R 



3ir , 
4 ' 



2R ^ 



. Choosing r such that 

(we may do so using the intermediate value theorem and the fact that R < r < h and h > 2R), 
we may determine that 




3tt / 5 / R R\ AR 
T 1 + 2\2h + 2hj + ~h' 



This implies 



(h) 



h) -)• 



R 



E(xl 



or 



E v 



.(h) 

'r 



E xl 



3vr 
" T 

37T /t 



as -| — > oo given that r = h. 



(12) 



(13) 



Remark 2. Recall that L = y |^4. j / vr anJ observe that Pr{h < a) < j^r. Therefore, we can 
obtain that Pr(h < u(N)) — > for a(N) = o(L) as N — > oo, which in return implies that 
Pr(h/R — y oo | rjdN + logrf — > oo as N — >■ oo) = 1. Hence, assuming that the conditions in 
TheoremUji hold, we have h/R — > oo a.s. as N — > oo. 

Remark 3. The asymptotic expected length of the routing path established by the random 1/2- 
disk routing scheme is = R/E (x' n ) ~ 2.36 times longer compared to the length of the routing 
path generated by the ideal direct-line routing scheme; in the ideal direct-line routing scheme 
we assume that there are relays located on the line connecting the source and destination with 
the maximal separation R, 

By averaging over all the possible source-destination pair distances h, we can determine the 
expected length of a typical random 1/2-disk routing path. Again using Pr [h < aR) < pr (aR) 2 
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and (fT2~l) we have that 



E (u R ) =E(E (vf | h) l {h < aR} + E (uf | h 



L{h>aR} 



< 



ira 3 R 2 



A 



T 1 



1 

to 2a 



and 



EM 



E E i/ 



ft 



+ 



37rE (ftljftxjR}) 



> 



R 



3vr /E(/i) 



/ 8a 



T* 1+4 



4 V i2 

The problem of quantifying E(/i) is well studied in the literature ifTSl . with the following 
known results for two network regions specifically: If the region is a planar disc with diameter 
2L, we have E (h) = 128L/(457r) w 0.9054L; and if it is a square with side length L, we have 
E(h) = (2 + V2 + 51og(V2+l))L/15 w 0.5214L. Choosing a = o (rf~ 1/6 ) and recalling 
Remark |2l we observe that E (hl{ h>aR y^ — > E (/i) as — > oo; hence, we have 



R , ^ 32 1 



(14) 



as — > oo. 



C. TheoremU\iii Proof: Variance of the Random 1/2-Disk Routing Path Length 

So far we have characterized the expected length of the routing paths generated by the random 
1/2-disk routing scheme. However, the expected value alone is not descriptive enough regarding 
the individual realizations of the routing path length. We need to determine how much the 
individual realization deviates from the expected value. Therefore, in this section, we consider 
the variance of the path lengths generated by the random 1/2-disk routing scheme. We first show 
that the variance is finite almost surely and then we show that asymptotically it grows linearly 
with the expected path length: 



Var(4 fe) |/Q VarQQ _ Qtt 2 



as iV — > oo. 

Consider the i.i.d. sequence {(x' n , y' n )} as defined in Section HV-Bl and define the generalized 
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stopping time Va to be : = inf{n : r n < a, r = b} for R < a < b < h. Observe that 
{ui b) > N} and {<} n< Ar are independent, and E (v^ j < oo and E ((x' n ) 2 ) < oo as shown in 
Section IIV-BI and Appendix El 
Note first that, by definition, 



CO A 



h- R< Yl (~9( r i-u x[, y'i)) = r - r u ( h ) An < h 



for any n, where An := min{z/j^\n}. Define 

U n := Xir=i( — 9(Ri x 'ii Hi))- F rom Wald's equation, Eq. ©, and the fact that g(r,x',y') is a 

nonincreasing function over r > R, we have 



E U 



CO . 
'An 

r CO a 
v R An 



h 



- E \Y1 (-9(n-i,x'i,y'i)) 



h \ < h 



i=l 



for all n. As shown in the previous section, it follows that E ( v R 



■iii 



h ) = linin^oo E ( Ujp A n 



h) < 



? < oo. Similarly, 

(E(-.g(i?,^,yO)) 2 Var(^ l) An | h) 



< 2 

< 2 

< 2 



Var 



(U w 

V V R 



Var {/ m 

v K „ 'An 



4/i i? 2 



+ Var {{yf A n)E (-.g(i?, ^ - A „ 
h) Var(-g(i?,^,yO) 



i? 4 



for all n, where the second inequality is due to Wald's identity ( ||20l , page 398). Thus, 



Var v 



Xh) 
R 



h 



lim Var ( u]t' A n 



^ 32h(h + R) 
' " R 2 <0 ° 



(16) 



This proves that the variance of path length generated by the random 1/2-disk routing scheme 
is finite almost surely. Next we will find some asymptotically tight bounds on the variance of 
the routing path lengths. We will frequently use the following well known inequalities 

VEp^)-v^In <y/E((X-Y)*), 
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and 



v/Var(X 2 ) - v/Var(F 2 ) < ^Var (X-Y). 



Let S v := Yln=i x 'n f° r a stopping time v such that {y > N} and {a4} n <iv are independent 
and E (u) < oo. Then by Wald's identity ($M, page 398) we have E (S v ) = E «) E (z/) and 

Var (uE «) — S w ) — E ((£„ - z/E (x'J) 2 ) = E (,/) Var «) . 

As such, we have 



'Var (z/)E (x' n ) - a/E (z/) Var (x' n ) 



In particular, for v = v^, we have 



'Var(i/E «)) - ^Var^E «) - < ^/Var^). 



\ 



Var 



1 4 h) i /i 



Var«) 
(E«)) 2 



< 



Var 


(s^ 1 /V) 




E 






(E (x 


J) 2 



(17) 



Hence, in order to prove the limit in (1151) . we need to show that 



Var ( S (h) I h 



E[u^\h) (E«)) 2 
as N —7- 00. Suppose R < a < b < h and note that 



Var S* w \ h 
. ^0 



lb 



-^(r^i,^,^) < x' n < -g(r n _ 1} x' n ,y' n ) + 



R 2 



2r n ,_i 



then together with (I8al) . we obtain 



v"'» 



S lb) — S (a 
V R v R 



U R 



^ X (-^K-l'^n.J/n)) + X 



R 2 



R 2 



,(*>) 



<b-a + R+ —v y R , 
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where the last inequality is due to the fact that r n > a for < n < v®. Therefore, we obtain 



Var (S w - S u 
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a, b 



S (b) — S (a) — E ( S (b) — S (a 
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< WE 
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V R V R 



a, b 
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2a R 
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< fl+ | E (^|6) + ^/Var(, s 
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< i? + 


2Kb 
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a 




5Rb 
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2' i ») 


< f a 



2a 



^326(6 + R) 



4" 



I < 32fc g+ fl ) . Finally, we let 



a, = i?(-|)^ fc , for = [log |f| and z = 0, 1, 2, • • • , k, where [log -|] is the smallest integer 
larger than log -| . Then we have 



From this, it follows that 



\ 



Var (s„m | fc) < J2 l/var (s„k, - 

fe 



i=l 



Oi-1 



{h\ lfk 

= 6kR + 5kRl-j 

< (6 + 5e)(l + log^) J R. 



Var f S v w | h 

R 

Rh 



h 



<(6 + 5e)(l + log^w|^0 
R \ h 



(18) 



as N — > oo, which concludes the proof for the limit in Eq. ([T5 
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Fig. 8. Random 1/2-disk routing realizations for A = 10 6 , \A\ = 1, and R = 5.2 x 10 3 , when the source is located at 
(1/3, 1/3) and its destination is located at (2/3, 2/3). 



V. Simulation Results 

In this section we compare our analytical results with some empirical results. In particular, 
Fig. [8] depicts some realizations for the routing paths generated by the random 1/2-disk routing 
scheme for an arbitrary source located at (1/3,1/3) and its destination at (2/3,2/3) with the 
following network specifications: \A\ = 1, A = 10 6 , and R = y^f^ ~ 5.2 x 1(T 3 . As illustrated 
in this figure, the path realizations do not closely follow the direct line connecting the source- 
destination nodes. The lengthes of the routing paths are 208, 208, 225, 223 for the realizations 
depicted in Fig. [8] (a), (b), (c), and (d), respectively. Based on (fT2l) we obtain the lower and upper 
bounds of 208, 256 for the expected path length with the asymptotic value of 211.3. (Note that 
the bounds derived in (|T2|) are for the expected path length; therefore, individual realizations for 
the path length might violate these bounds.) 

In Fig. |9l we compare the (normalized) empirical mean, -^E^z/^j, of the path lengths 
generated by the random 1/2-disk routing scheme with the analytical bounds derived in Eq. 
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Fig. 9. Numerical comparison between analytical bounds derived in Eq. J 1 2b and the (normalize d) emp irical mean of the path 
length generated by the random 1/2-disk routing scheme when h — \/2/3, \A\ = 1, and R = -J °f ■ 



(fT2l) . As shown in this figure, the normalized empirical mean converges to 3n/A pa 2.3562, and 
is always bounded by the expressions derived in Eq. (fT2]> . 

In Fig. [TOl we compare the em pirical standard deviation for the path length normalized by 
the mean, y^Var (^r^ /E (vr^, with the analytical bounds derived in Eq. (fT8l) when the 
source and destination are h = a/2/3 distance apart and the transmission ranges are chosen as 
R = y // 21 ° gA for different values of network node density. As shown in this figure, the normalized 
empirical standard deviation converges to a/97t 2 /64 — 1 pa 0.6228, and is always bounded by 
the expressions derived in Eq. (fT8l . Furthermore, it can be seen that although the bounds in (fT8l 
are quite loose for small values of A, the asymptotic standard deviation derived in (fT5l) is very 
close to the empirical standard deviation even for small values of A. 

In Fig. [HI we demonstrate the deviation of the path length realizations from its asymptotic 
expected value when the source and destination are h = y/2/3 distance apart and the transmission 
ranges are chosen as R = .y/ 21 ° gA for different values of network node density. As shown in 
this figure, the deviation of the path length realizations increases as the network density and 
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Fig. 10. Numerical comparison between analytical bounds derived in Eq. J 1 7b and the (normalized) empirical s tandard deviation 
of the path length generated by the random 1/2-disk routing scheme, when h = \/2/3, \A\ = 1, and R = 2 '° s A . 

consequently the expected length of the routing path increases. However, all realizations stay 
relatively close to the value predicted by Eq. (fl"3l) . 

VI. Generalization 

In the previous sections we derived sufficient conditions for the network to be connected de- 
ploying the random 1/2-disk routing scheme and quantified the mean and variance asymptotes of 
the routing path generated the random 1/2-disk routing scheme. In this section we present some 
guidelines that generalize the aforementioned results for some other variants of the geometric 
routing schemes such as MFR, DIR, NFP, and the random 77-disk routing scheme, where the 
latter one is the generalized version of the random 1/2-disk routing scheme with an ?7-disk as 
its RSR. 

Observe that the results of Section [III] were derived for the general ??-disks relay selection 
region which encompasses most of the geometric routing schemes such as MFR, DIR, NFR 
and the random //-disk routing scheme. Let A be the set of all nodes (in the RSR of a specific 
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Fig. 11. Random 1/2-disk routing realizations for A = 10 6 , |A| = 1, and R — 5.2 x 10 -3 , when the source is located at 
(1/3, 1/3) and its destination is located at (2/3, 2/3). 

transmitting node) that can be selected as the next relay by the relay selection rule (RSL) of 
the geometric routing scheme. For example, in the cases of MFR, DIR, NFR and the random 
77-disk routing scheme we have: A M fr : = {(x' n ,y' n ) £ |RSR : x' n > x, for all (x,y) G |RSR}, 
A D ir := G |RSR : | tan' 1 (y'J x' n )\ < \ tan-\y/x)\, for all (x,y) G §RSR}, A nfp := 

{(x' n ,y' n ) G |RSR : ^(O 2 + (Vn) 2 < \/(^) 2 + (v) 2 , for all (x,y) G |RSR}, and A, = 
{(x' n ,y' n ) G r/RSR}, respectively. Since the nodes in A (if more than one) are indistinguishable 
by the RSL, the transmitting node selects one of the nodes in A randomly as the next relay. 
Next, we present the generalized results for the network connectivity and the mean and variance 
asymptotes of routing paths generated by the general geometric routing schemes. 

Corollary 1. Let A be the set of all nodes that can be selected by the relay selection rule as 
the next relay. Then the network is connected employing the geometric routing scheme a.a.s. if 
E(g(R,x',y')l {A} ) < 0. 
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Proof: The proof is immediate due to ( TTOT ). ■ 

Corollary 2. If E [g(R,x',y')l{ A }) < and E ({y') 2 l{A}) < RE {x'1{a}), the expected length 
of the routing path generated by the general geometric routing scheme connecting a source- 
destination pair that is h-distance apart scales as E{y | h) ~ h/E (x'l {A} ) as N — > oo. 

Proof: The proof follows directly from CCD) and noting that if E ((?/) 2 l{A}) < -RE (x'l{ A }), 

using the intermediate value theorem, we can find r such that / 2 ' R ^~ fe - > \ = ? , R — r (\ Urn + 1), 

E{(y'yi{A}) E(z'l {A} J v 2H 

which yields the bound in Eq. (fT2l) and hence the desired result. ■ 

Corollary 3. IfE (g(R, x', y')l{A}) < 0, the variance of the path length generated by the general 
geometric routing scheme, normalized by its mean, scales as Var{y) /E (z/) ~ Var (x'l^}) I (E {x'l 
as N oo. 

Proof: The proof follows the same steps as in Section IIV-C1 ■ 



VII. Conclusion 

In this paper, we presented a simple methodology employing statistical analysis and stochastic 
geometry to study geometric routing schemes in wireless ad-hoc networks, and in particular, 
analyzed the network layer performance of one such scheme named the random 1/2-disk routing 
scheme. We defined a notion of network connectivity considering the special local properties 
of geometric routing schemes and determined some sufficient conditions that guarantee network 
connectivity when each node finds its next relay in the so-defined 1/2-disk. More specifically, if 
all nodes transmit at a power that covers a normalized area d and the expected number of nodes 
in the network is N, the network is connected a.a.s. if d — o (iV~ 2//3 ) and rjdN + logd — > oo 
when N — > oo. Furthermore, we showed that the process of path establishment by the random 
1/2-disk routing scheme can be approximately characterized by a Markov process that converges 
statistically to the actual process asymptotically. Then using this Markov characterization, we 
derived exact asymptotic expressions for the mean and variance of the path length. Furthermore, 
we provided guidelines to extend these results to other variants of geometric routing schemes 
such as MFR, DIR, and NFP. 
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Appendix A 
Proof of Proposition [j] 

First, let us consider the distribution of a Poisson point process conditioned on deleting one 
point. Let $ be a homogeneous Poisson point process with intensity A and assume a fixed region 
D. If 4>{D) > 0, one point in D is selected at random and removed. Let X be the location of 
that point. The distribution of $ on D c remains Poisson and independent of $ on D, and thus 
independent of both 4>(D) and X. Let $' be the (point) process with the point at X deleted. 
(Note that the distribution of $' is not the same as the reduced Palm distribution j[2T| of <3>, as 
the location of node X is random.) 

Let Ai,A 2 ,-- - , A k be a partition of D. Given 0(D) > 0, the points in D are distributed 
uniformly. If one point is removed at random, the remaining points are still distributed uniformly 
on D. Hence, 

* ((VOW = n j} | m > 0, X )= (1 - e-n- £ ft 

\j=l / i=l j=l J \J >' 

I A \\rij 

-\\A 3 \ 



\\D\ Ami) 



e -x\D\)( ni + --- + n k + l) n ' ( nj ")l e 



(19) 



since |Ax| + • • • + \A k \ = \D\. Therefore, conditional on 4>(D) > 0, $' is independent of the 
location of the removed point (X). In particular, 

= Pr (j>(D) = n + 1 0(D) > o) . 
Furthermore, given ni + ■ ■ • + n k = n > 0, we have 

\ k / I A I 

|£>l 



Pr ^Q{0'(A) = | > =n,Xj = _ n _ n J n 

Thus, for A C D and given (j)'(D) = n > 0, is conditionally Binomial ^n, j^jj- Without 
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knowing (j)'(D), however, we obtain from (fT9l) that 

/ \ AlDle _A|z?l 

Pr(f(A) = k <J>(D)>0,x)=-± 



E 



(X\A\) k (X\A C n D\y 



e -x\D\\ ^ k + j + 1 k\ 

1 3=0 J 



Al^le-^l (X\A\) 



k /•! 



'l- e -MDV 



k\ 



yk e W™>\V d y 



(20) 



where the second equality is due to 



f 1 a> = V— [°—ix= [" ^-e'ix = f /e-cto 

After the aforementioned preliminaries, we now proceed with the proof of Proposition [Q 
Suppose C is a random set that depends only on The points of if any, which are in 
CD := C D D, are uniformly distributed and independent of the points in CD C , which are also 
uniformly distributed (if any). The combined points are uniformly distributed on C only if the 
expected proportion of points in CD is ^3-. 

However, the expected proportion of points in CD is strictly less than in our case as we 
now compute. Given 4>'{C) > 0, the probability that a randomly selected point in C is also in D 
is E | <t>'{C) > O,0(L>) >0,x\ Let = when (//(C) = 0. Using ©, we have 

Pr(V(C)>0 <l>(D)>0,x\=l-PT(4>'(CD) = 0,<l>'(CD e ) = Q </>(D)>0,x\ 

\\CD C \ 



1 - 



X\D\e 


-A|D| g A|C c D| _ 1 


1 - e- 


A|C C L>| 6 


|D| 


1 e r -A|C| 


|C C L>| 


l_ e -A|D| 



so we have 



\D\ 



\C C D\ 



A l c l < Pr U'(C) > <j>(D) > 0, x\ < 1 - e- A l c l 



Using the observation above and (1201) we obtain 

* 7 n— m— 1 u 

l) dw J 1 y n e^ CCD ^dy 



\\ D \ e -x\cuD\ ~ (A|OD|) n r 1 n-1 



l_ e -A|D| 



^ (A|C£>|) n r 



w e 



,A|CZ5 c |iii 



13 Note that D and C here correspond to D n -i and D n in Section ITVl respectively. 
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\\D\t 



-X\CUD\ 



1 - e-^l 



\\CD\ye^ CD \y w+ \ CCD ^ (e^ CD ^ ~ l) AwAy 

(21a) 



< A|g|e ^ I' X\CD\ye^ CD \y w +\ CCD \y+\ CDC ^d W dy 



l- e -x\D\ 
X\D\ e -^ CuD \ f 1 \CD\y 



< 



l_ e -\\D\ J Q \CD\y+\CD 

\CD\ X\D\e- x \ CuD \ f 1 
\C\ 1-e- 



AID 



^ e X(\CD\y+\CD"\) _ A e X\C c D\y iy 



f e \QCD\„+\CD°\) _ A e \\C°D\y dy 



\CD\( \D\ , C | 

\C\ \ \C C D\ 

\CD\ 
\C\ 



-Pr 



(V(C) > > 0,x) . (21b) 



Therefore, 



Noting that 



0'(C) > 0, <£(£>) > 0,X) < 



\CD\ 
\C\ 



1 - - < ae" a / ?/e as/ dy = 1 
a 



1-e" 



< 1 



we could derive (l22l) from (I21al) 

W ' (c)>0} I ^ > °' J = i _ e -A|£>| 

- / / \\CD\ye* 
Jo Jo 



\CD\y 



o \CD\y + \CDc 



,X{\CD\y+\CD°\) 



1 e 



A|C c _D|y 



dy 



\CD\yw + \C c D\y) 



dwdy 



> 



\\D\e 



-A|cun| 



l- e -x\D\ 
\CD\ 1 



1 \CD\ 
\C\ ■ 



(eH\CD\y+\CD°\) _ ^ e MC<D\y dy _ f e X\C°D\y ^X\CD\y _ &y 



> 



> 



> 



\c\ 


(1 


- e -\\D\) 


\c\ 


( 


_ e -A|CD«| 


\CD\ 






\CD\ 




1 


\C\ 


(1 


_ e -A|r>|) 


\CD\ 




1 


\C\ 


(1 


_ e -A|£>|) 


\CD\ 




1 


\C\ 


(1 


- e -MD\) 



\\D\e- x \°\ / ye x Wydy-\\D\e- x \ CuD \ 



II'' 



\C\ 
\CD\ 



-\\CD C 



1^1 

\C C D\ 



£ -X\C\ _ e -\\CUD\ 



1 



|C| 



X\D\ \CD\ 
1 



-\\CD a 



\C C D\' 



-A|cun| 



\C\\D\ e - 2 x\CD° 



(22) 



for large enough N such that 1 — — ^wj^pj exp(— 2A|CT> C |) > 0. Hence we can ascertain 
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that 



U\CD) \ f 1 \C\\D\ 2Xlcm \ \CD\ 

As such, the selected point is less likely to be in D than the case where we assume $' is 
Poisson on C. 

Appendix B 
Derivation of inequality ® 

We have «, y' n ) = (Rvcos(9), Rvsm(6)), where 6 ~ Uniform (-vr/ 2, vr/2) and v ~ Beta(2, 1) 
are independent. Thus, we have 

2 f n / 2 AT? 

B(x ; n ) = R- cos(9)d6 / 2v 2 dv = — , (23a) 
7t Jo Jo 3tt 

Z{(y' r f) =— r /2 sin\6)d6 !\iMv = ^-. (23b) 

n J-n/2 JO 4 

Also, by first changing x to 1 — x and then using polar coordinates, we obtain 

^E(g(R,x' n ,y' n )\ +1 = ^ J l x 2 +y i<i\f {I - x) 2 + y 2 dxdy 

= - ( [ l(i_ :E )2 + j / 2< lV / x 2 + y 2 dxdy 
71 Jo Jo 

c\ pit/ 4 rsccd 

= - I I 2v 2 dvd9 + 

7T 



JO 



2 rn/2 /.2cos6» 

2v 2 dvd6 



71 



tt/4 JO 



4 r^/ 4 
3tt 7o 



secfl) 3 + (2sin#) 3 ) dfl 
3(2 3 / 2 ) + 6 log(l + v^2)) + 64 - 5(2 7 / 2 ) 



9tt 

0.7499728 . 



Hence, E (g(R, x' n , y' n )) < - 
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